Theory of Graphene Chiral Quasiparticle LDOS maps 
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We present a theory of momentum space local density-of-states (LDOS) maps N(q, uj) in graphene. 
The LDOS map has both intravalley contributions centered near zero momentum and reciprocal 
lattice vectors and intervalley contributions displaced by the wavevector K' — K which connects 
graphene's two distinct Dirac points. Using graphene's Dirac equation chiral quasiparticle continuum 
model, we obtain analytic results which explain the qualitative differences between these two LDOS 
map features. We comment on the sensitivity of both N(q,cj) features to the mix of atomic length 
scale and smooth disorder sources present in a particular graphene sample. 



I. INTRODUCTION 

Graphene, a honeycomb lattice of carbon atoms, is a 
recently realized^ two dimensional (2D) electron system 
(2DES) with a variety of unique properties^ This in- 
triguing electronic system is now being actively explored 
both experimentally and theoretically. Graphene is de- 
scribed at low energies by a 2D masslcss-Dirac wave 
equation in which the role of spin is assummed by a 
pseudospin which represents the two atoms in its unit 
cell. Many of the unusual properties of graphene sheets, 
including the shift in the densities at which quantum 
Hall^i pleateaus occur, follow from momentum-space 
Berry phases^ associated the pseudospin (sublattice) 
degree-of-freedom. Graphene quasiparticles have defi- 
nite pseudospin chirality, i.e. definite projection of pseu- 
dospin along momentum measured from one of the two 
independent Brillouin-zone corner points at which the 
gap vanishes. The experimental observation of the quan- 
tum Hall effect in graphene was an important demonstra- 
tion that these 2DESs behave, at least in some respects, 
very nearly in the ideal way anticipated by Wallace^ 
many years ago in his early analysis of the electronic 
structure of graphite. Wallace realized that the states 
near the Fermi level of a graphene sheets should separate 
carbon 7r-orbital bonding and antibonding bands, and 
that the gap between these bands would vanish at two- 
points in the honeycomb lattice Brillouin-zone, the Dirac 
points of graphene. Experiments have confirmed this 
simple picture in most respects, something that was not 
a priori obvious given the potential for either electron- 
electron interactions or disorder to alter physical prop- 
erties. Indeed the quantum Hall effect, which can be 
viewed as a topological property of the 2DES, tends to 
be forgiving on details. In view of the tremendous inter- 
est in studying graphene sheets and characterizing their 
disorder, there is strong motivation for expanding the 
comparison between experiment and theoretical analysis 
based on Wallace's 7r-band model to new observables. 

In this article we present theoretical predictions for 
the local density-of-states of weakly disordered graphene 
sheets, a quantity which can be measured using scanning- 
tunncling-microscopy (STM). A new appreciation of the 
ability of this type of measurement to shed light on the 



character of the disorder in a sample, and also on un- 
derlying clean system electronic properties, has emerged 
from a highly successful series of studies of cuprate 
superconductors^. The experiments rely on the ability to 
make stable atomic-resolution STM scans of the LDOS as 
a function of energy to over a large real-space field of view. 
(The energy oj is varied by changing the bias voltage be- 
tween the STM tip and the sample.) Preliminary exper- 
imental data on graphene sheets is already available^, 
and we can expect that further refinements in sample 
quality and experimental technique will enable detailed 
analysis which will extract much useful information. Pre- 
vious theoretical wor k 11 ' 12 ] 13 ! 14 ? 15 has discussed the nu- 
merical construction of LDOS maps from graphene's hon- 
eycomb lattice 7r-band tight-binding model. Analytic ex- 
pressions for the LDOS modulations were obtained by 
Bena*i2 through an expansion in powers of 1/r of the 
amplitudes in real space. Peres et al. have looked at 
impurity induced localized states at low energies in the 
T-matrix approximation™ In this theoretical contribu- 
tion we focus on the Dirac equation continuum limit of 
the 7r-band model, from which it is possible to obtain an- 
alytic results which we believe can contribute to a more 
meaningful interpretation of experimental LDOS maps. 

In this paper we suggest that it can be useful to mea- 
sure A and B sublattice LDOS maps separately, some- 
thing which is possible in principle since the experiments 
have atomic resolution. As we explain, the difference sig- 
nal (which will normally be weak) and the sum signal 
are complementary probes of a sample and its disorder. 
With this method one is able to extract more informa- 
tion and consequently identify the type of disorder in the 
system. 

Because of the disorder always present in real mate- 
rials the LDOS never has the lattice periodicity. When 
Fourier transformed to momentum space, the additional 
modulations^ reveal the dominant momentum transfers 
between quasiparticle states at energy uj (measured from 
the Dirac point) and hence map out the electronic struc- 
ture. In a superconductor the Bogoliubov-deGennes co- 
herence factors, which specify the particle and hole am- 
plitudes in BCS quasiparticles, help determine the rela- 
tive probability of scattering between different constant 
energy surface segment al 19 . As we shall see shortly, 
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the pseudospin mixing present in graphene quasiparti- 
cle excitations plays a similar role. At small uj the 
momentum-transfers in graphene arc naturally classified 
as intra-valley (small momentum transfer close to the 
same Brillouin-zone corner) or inter-valley (large mo- 
mentum transfer from one valley to the other.) We will 
see that both intra-valley and inter-valley scattering am- 
plitudes are influenced by pseudospin chirality, but be- 
cause the two valleys have opposite chirality their STM 
momentum-space maps differ qualitatively. 

Our paper is organized as follows. In Section II we 
briefly summarize the weak-scattering analysis which has 
been used to provide STM momentum-space maps with 
a simple and consequential interpretation. In Sections III 
and IV we apply the analysis first to a graphene 7r-orbital 
tight-binding model and then to the low energy Dirac 
equation continuum model. The Dirac equation model 
allows many elements of the calculation to be carried 
out analytically, enabling us to provide more guidance on 
the qualitative interpretation of STM momentum space 
maps. We conclude in Section V with a brief summary 
and discussion. 



II. THE WEAK DISORDER APPROXIMATION 

A. Real-space LDOS maps 

The influence of a scattering potential V(r) on the 
LDOS of an otherwise clean system can be described by 
expanding the Greens' function in powers of V(r). LDOS 
maps are most revealing when the scattering potential, 
V(r), is weak, justifying truncation at first order - the 
Born approximation. For a simple parabolic band, this 
approximate treatment leads to the well known Freidel 
oscillations^. In a LDOS-map experiment, weak disor- 
der provides an electronic system with its own weakly- 
coupled probe - one which is able to provide momentum- 
resolution of single-particle properties. It is generally as- 
sumed that the weak-coupling approximation is at least 
qualitatively valid whenever sharp features appear ex- 
perimentally in the momentum-space maps we describe 
below. Maps similar to the ones we calculate later in 
this paper should emerge from LDOS-map studies of 
graphene if samples of sufficiently high quality can be 
prepared. 

To model graphene we consider the following per- 
turbed Hamiltonian: 

T~t = Tto + H. v 
H v = Y,V(r)n(r) (1) 

r 

where the unperturbed Hamiltonian TIq is either the it- 
band tight binding 8 -^ model for graphene or its Dirac 
continuum limits The Dyson equation for the Greens' 



function is then 

G(r,r',oj) = G"(r'-r,uj) 

+ J dsG°(s- r, lu)V(s)G(s, r', w), (2) 

where G°(r' — r,ui) is the unperturbed Greens' func- 
tion of a clean-limit graphene quasiparticle at energy hw 
and G(r,r',oj) is the inhomogeneous disordered system 
Greens' function which includes the effect of scattering 
by the potential V(r). In order to describe a non-Bravais 
lattice, it is convenient to write the Hamiltonian and 
therefore the Greens function as 2 x 2 matrices with the 
argument r defined only on the Bravais lattice. We de- 
fine the vector r to lie on the A sublattice and the B 
sublattice atoms are obtain through a shift by the vector 
r. For on-site disorder V is diagonal and its elements 
specify the A and B site potentials in a particular unit 
cell. It is convenient to write 

V(r) = Vo(r)a + V 3 (r)a 3 , (3) 

(Here V + V 3 = V A , V - V3 = Vb , (Jo is the 2 x 2 identity 
matrix and 03 is the diagonal Pauli matrix.) Disorder 
sources which are smooth on an atomic scale, for ex- 
ample smooth ripple a 24 i 25 ' 26 or Coulomb potentials from 
remote ionized impuritie s 27 ' 28 ! 29 , will contribute only to 
Vo, whereas atomic scale disorder sources^ like impurity 
atoms or vacancies, will have large differences between 
Va and Vb and contribute to both Vb and V3. The local 
density of states at position r and energy hu, N(r,u>), is 
given by 

= _l x ilm[G AA (r,r,u J )] on A site 
7r 1 Im[Gflfl(r,r,w)] on B site 

B. Momentum-space LDOS maps 

The zeroth order term in the potential expansion of 
the Greens' function gives the clean-limit LDOS which 
is periodic and therefore has non-zero Fourier compo- 
nents only at reciprocal lattice vectors. The disorder- 
induced non-periodic spatial modulations appear in the 
subsequent terms. The Born approximation sums over all 
quantum paths which include a single scattering event. 
In terms of the exact single-particle eigenstates of a disor- 
dered graphene system, the Born approximation to the 
LDOS accounts for the corrections to the Bloch wave- 
functions of the perfect crystal which appear at first order 
when disorder is treated as a perturbation. In principle 
one can continue the perturbation scries to include higher 
orders of the potential. For a single-scatterer, the full se- 
ries is easily summed to construct the full T-matri x 12 ' 31 , 
a necessity when strong scattering leads to impurity res- 
onance states. In this paper, motivated in part by the 
absence of any experimental evidence for resonances in 
typical graphene sheets and also by the possibility of ex- 
tracting more information about the sample from STM 
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experiments when this limit applies, we concentrate on 
the weak disorder limit. 

A STM measurement naturally projects a quantum 
wavcfunction to one sublatticc or the other and not, 
as implicitly assumed by Bena and Kivelsonii, to a k- 
dependent Bloch pseudospinor. This point is discussed 
in detail by Bena and Montambaux^. We restrict our 
attention, without loss of generality, to evaluating the 
LDOS measured on the A sublattice. This sublatticc 
separation will allow distinction between the two types 
of scattering profiles as will be demonstrated shortly. In 
the Born approximation, we may restrict our attention 
to a single-site impurity potential which we take to be 
located in the unit cell with lattice vector L = 0; the 
distributed potential case can be constructed simply by 
adding contributions from different unit cells. The first 
order correction in the Greens' function is^ 

G« (r, r, w) = [G°(r - 0, u;)V(0)a 0/3 G (0 - r, cu)] AA = 

V(0)[G° AA (r)G%(-r)±G AB (r)G° BA (-r)}.(5) 

In Eq.(f5|) the first term describes the case of measuring 
and scattering on the same sublattice, and the second 
term the case of measuring on one sublattice and scat- 
tering from the other. 

Electronic structure information is revealed most di- 
rectly by Fourier transforms of the LDOS map. We 
evaluate the momentum-space LDOS map first from 
graphenc's 7r-orbital tight-binding model and later take 
the continuum limit of these calculations to obtain ana- 
lytic Dirac- model results. The tight-binding model Bloch 
Hamiltonian 

Uak = U 7 o ) 

Ik = t(e lk y/^- + 2e lfc «/ 2v/ 3cos(fc x /2)), (6) 

where t is the hopping amplitude and wc measure length 
in units of the triangular lattice constant (not the carbon- 
carbon distance) a=2.46A The 2x2 unperturbed Mat- 
subara Greens' function is then (iui — 7i) _1 . 

The tight-binding model Hamiltonian (Eq.©), and 
hence G , is constructed in a representation in which 
Bloch state basis functions have a phase difference k ■ r 
between site 2 and site 1 in every unit cell. Here r is the 
vector from site A to site B. It is critical that the same 
gauge is used when the Fourier transform of the LDOS is 
constructed, in particular in evaluating the Fourier trans- 
form of the LDOS on one sublattice due to scattering 
from a site on the other: 

J2e- i ' lR G AB (R,u J )G BA (-R,w) = 

R 

Y,e~ lRq e lk(R+T) GAB(k,u)e-* k ' (R+T) G BA (k\u;) 

R fc,fc'GBZ 

= £ e^-^G%(k^)G BA {k-q + Q,u){7) 
feeBZ 



The exp(i(fe — k 1 ) ■ r) in the second line results from the 
gauge choice in which the origin is taken to be on the A 
sublatticc. In the last line wc have noted that the lattice 
vector sum yields a 5-function which restricts k — k' to 
be equal to q up to a reciprocal lattice vector Q. Q is 
therefore chosen such that k^q + Q is in the momentum- 
space primitive cell. The full Fourier transformed LDOS 
is therefore given by: 

5N A (q,u) = ^-Y,V^q)(A fl (q,u-iS)-A^q 1 u + iS)) 
M=0,3 

A^g.iw) = Y.\- G0 AA(k,iLo)G AA (k-q,iu) (8) 
k 

± e iQr G AB (k,iu)G° BA (k-q + Q,iu)} 

where the retarded and advanced parts of A are obtained 
by analytically continuing ito — * to ± i5 and V^{q) is the 
Fourier transformed potential. Note that the factor e~ tqT 
which appears in Eq.([7]) is not present in Eq.©. This 
is due to a 'form factor' which is implicit in the Fourier 
transform of the potential V(q), i.e., due to the distance 
r between the two atoms in the unit cell a phase of e iqn ~ 
appears in the B-componcnt of the potential. These ex- 
pressions do not include the 7r-orbital form factor which 
is expected to gradually decrease momentum-space am- 
plitudes at wavevectors outside the Brillouin zone. 

It is possible to describe any arrangement of impuri- 
ties from Eq.Q by combining the results for uq and 03 
potentials. For example, a single impurity on the A sub- 
lattice is represented by ctq + 03 (and therefore setting 
Vb = V3) . It is also possible to obtain the amplitude of 
the LDOS modulations on the B sublattice from Eq.© 
by replacing labels A <-> B and reversing the direction of 
the basis vector, r — > — r. The single impurity results 
are given by: 















-A#( 9 )) 




Af(<7) = i(A B + Af) 


= la^r 


- (A3O*) 




Al(<7) = i(A B -Af) 


= la^r 


+ (A#)*) 


(9) 



where the subscript represents the position of the impu- 
rity and the superscript is the sublattice on which the 
measurement is done. 

In this work we have assumed that the experimen- 
tal data is separated to the two sublattices before the 
Fourier transform is performed. If the experimental data 
is Fourier transformed without the sublattice separation 
the two diagonal components of the perturbed Greens' 
function should be added (G AA + Gbb)- In the case of a 
single impurity (as studied by Bena—) the result obtained 
in this way is identical to our 5N A (q). In the case of two 
identical impurities on the A and B sites, the FT con- 
structed from both the A and B sites (SN A+B ) is equal 
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to the real part of our SN^. In the case of two equal and 
opposite impurities (173 scattering) the combination of A 
and B sublattices cancels the real part and the resulting 
5N^ +B (q) is purely imaginary and is equal to the imag- 
inary part of our 5N^-(q). These identifications of real 
and imaginary parts apply, of course, when only one unit 
cell of the lattice has scatterers, and the A site of this 
sublattice is chosen as the origin of coordinates. In the 
more general case of distributed disorder, which we ex- 
pect applies to real materials, it is nevertheless true that 
the difference between A and B sublattice maps arises 
purely from the A3 response whereas the sum of the A 
and B maps will be dominated by the Ao response if the 
dominant disorder varies smoothly on an atomic scale. 



III. TIGHT-BINDING MODEL NUMERICAL 
RESULTS 



Real Imaginary 




FIG. 1: The real and imaginary parts of 8N / 3 (q) repre- 
sented in grey-scale (dark is high) in momentum-space. The 
LDOS maps were constructed from graphene's 7r-orbital tight- 
binding model. These numerical results were obtained at en- 
ergy ui = 0.2t. The solid lines (red online) are a guide to the 
eye connecting reciprocal lattice vectors. The dashed lines 
(blue online) are a guide to the eye connecting inter-valley 
features in the LDOS map. The dashed lines also form the 
triangular lattice Brillouin-zone boundary. 

Momentum-space LDOS maps constructed by evalu- 
ating Eq. © from the full tight binding Hamiltonian are 
illustrated in Fig.©. These results are for identical am- 
plitude scatterers, the SN (q) result, and opposite scat- 
terers, the dNs(q) results, in the unit cell at the origin. 
Results for a general scatterer can be obtained by in- 
serting its Fourier transform in Eq.©. Efficient evalua- 
tion of the LDOS is achieved by Fourier transforming the 



unperturbed Greens function to real space (through the 
fast Fourier transform algorithm), multiplying matrices 
appropriately, and then taking the imaginary part. This 
quantity is then Fourier transformed back to momentum 
space, mimicking the procedure used to analyze exper- 
imental data. A closer zoom on intra- and inter-valley 
features from Fig.© is given in the inset of Fig.© and 
in Fig.© respectively. 

These results arc filled with interesting features which 
reflect the physics of graphene sheets, motivating the ex- 
perimental studies which this work aims to assist. Most 
obvious is the expected appearance of clear separate fea- 
tures associated with intra-valley and inter-valley scat- 
tering. Bloch states near the Fermi energy of neutral 
graphene sheets appear close to the two valley points, 
K = (4tt/3,0) and K' = (8tt/3,0) at which the tt bond- 
ing and antibonding bands meet. States with an energy 
ui (measured from the neutral system Fermi level) oc- 
cur close to a circle centered on K or K' with radius 
k u = lu/v where v is graphene's Dirac-cone velocity. 
Smooth disorder potentials contribute only to Vo(q) and 
have large amplitudes only for scattering within these val- 
leys. They therefore contribute to LDOS map features 
only near q = or reciprocal lattice vectors. We see in 
Fig.© that the intra-band A (q) features at \q\ = 2fc w , 
associated with scattering across a Dirac cone, are much 
weaker than the corresponding features in As(q). We also 
note that the inter-valley features which appear in the 
LDOS map near wavevectors ±(K—K)' have an interest- 
ing angular variation which is absent in the intra-valley 
feature. Neither intra-band or inter-band features are 
periodic under translation by a reciprocal lattice vector, 
as explained previously. Another feature of the LDOS 
modulations which appears in the numerical results is 
the structure in the vicinity of reciprocal lattice vectors, 
i.e., in the corners of the (red online) solid hexagonal 
zone in Fig.©. The amplitude of the modulations near 
these points may be obtained from the modulations of 
small momentum transfer through a 2n/3 rotation of the 
B-atom scattering contribution. The above features may 
be understood in the context of the Dirac model, pre- 
sented in the next section. 



IV. DIRAC MODEL ANALYTIC RESULTS 

One of the reasons for the excitement around graphene 
is its low energy behavior. At two valley points in the 
Brillouin zone the energy bands touch and the dispersion 
is linear. This makes graphene a zero gap semiconduc- 
tor whose low energy Hamiltonian is a condensed matter 
realization of the Dirac model. Using the massless Dirac 
model, we can achieve a deeper understanding of the nu- 
merical results reported in Fig.©. 

For our analytic calculations it is convenient to choose 
a unit cell in momentum space which includes both val- 
leys in its interior. The more symmetric triangular lattice 
Brillouin zone (BZ) has the disadvantage of separating 
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FIG. 2: FT LDOS of A (dashed, blue online) and A 3 (solid, 
red online) type scattering, as a function of \q\ in unites of 
2uj/v, around zero momentum transfer. Note that the FT 
is purely real and dependent only on \q\. The inset shows 
the functions in momentum space (to be compared with the 
numerical data presented in Fig.fl]). 



each valley into three pieces which together appear at 
the six BZ corners. We choose instead as a unit cell the 
parallelogram constructed from the two reciprocal lattice 
vectors: a\ = (2n, 27r/v / 3) and a-i = (2n, — 2n/\/3). The 
valley points are then K = (47r/3, 0) and K' = (8ir/3, 0). 
Linearizing the tight-binding Hamiltonian around these 
points leads to: 



njDirac 

Ati. 





I ^ ^^1] 







(10) 



where the (+)- sign corresponds to the (K')K-valley and 
v = \/~3t/2 is the quasiparticle velocity. The difference in 
sign corresponds to a difference in sign in the pseudospin- 
chiralit y 34 ' of the Bloch states at the two- valleys which 
is responsible, as we explain below, for many of the qual- 
itative features of the LDOS maps. The unperturbed 
Matsubara Greens' function is therefore given by 



G°(k,ico) 



1 



k 2 



±k x -iK -v. )• ( n ) 



lul 



For convenience we have rescaled our momentum by the 
velocity v. We may now evaluate Ao and A3 within the 
linearized model. 



A. Intra-valley Scattering 

Here we use the same Greens' function for both the 
initial and final states (i.e. the same chirality sign) and 
arrive at: 



A / 3 (<7) = 



d 2 k —cj 2 ±(k x + ik y )(k x 



(2tt)= 



{to 2 + k 2 ){to 2 + (q - k) 2 ) 



(12) 



where the (-)+ sign corresponds to (As)Ao. The 
right hand side of Eq. (fT2|) can be simplified with the 
Schwingcr/Feynman trick^: 



iky)(k x 



iky + iq y ) 



1 f d 2 k -iv 2 ±(k 

3 J (2tt) 2 [k 2 + lj 2 + (1 - x)q 2 - 2k ■ q(l - x)) 

1 r d 2 k -u? ± (k 2 - x(l - x)q 2 ) 

d X J W? (k 2 + A) 2 



(13) 



where A = tu 2 + x(l — x)q 2 . Integrating over momenta, 
we find that the terms with momentum independent nu- 
merator yield (— lo 2 =p x(l — x)g 2 ))/47rA, whereas the k 2 
numerator terms yields 1 + log(A/iT))/47r where f2 is an 
ultraviolet cutoff that arises in the dimensional regular- 
ization scheme 3 ^. Integrating over x we find that for Ao 
the terms with momentum independent numerators give 
— 1/4-7T so that 



-5 



dx log 



1 



-1 



4tt 



2uj 



uj 2 + x(l — x)q 2 
2 



71 2f + CT + log - (14) 



where 



T[z) = \/— ' /2 — 1 arctan 



The physical quantity is given by: 



5N A (q) 



2^ 2 



Im 



V-z 2 - 1 



2iu> 



(15) 



(16) 



It is interesting to note that T vanishes when |q| = 2oj. 
This is despite the fact that energy conservation leads 
to the requirement that the initial and final states be 
on the same contour of constant energy with radius ui. 
This property requires the (-) sign in the fc-independent 
term to cancel a singularity at |q| = 2uj; the singularity 
does appear when the opposite sign is taken in evaluating 
A3. The absence of this singularity in the intra-band Ao 
map is largely responsible for the qualitative difference 
between intra-band and inter-band maps. 

The physics behind this cancelation can be understood 
qualitatively as follows. We may evaluate the dominant 
contribution to A /3 in terms of scattering between eigen 
states of the unperturbed system. These states are 2- 
vectors (pseudospinors) of the form ^ k = (l,±e~ l ^ h ) 
where the (-)+ sign corresponds to the (negative) pos- 
itive energy band and <pk = axctaxi(k y / k x ) near the K 
valley. Due to energy conservation and to the sharpness 
of quasiparticlcs, the dominant scattering events involve 
on-shell states with energy lo. The phase space integral 
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for elastic scattering on the energy shell. 



d 2 qS(uj - v\k\)S(uj - v\k - q\) 



Real 



Imaginary 



dqqdO S(uj - [(lj + vq cos(6>)) 2 
d(cos(0)) r/ ___ /m , q 



v 2 q 2 sm 2 (e)} 1 ' 2 ) 



yj\ -cos 2 (6>) 
2luQ(2lu- vq) 
v ^/Tu 



6(cos(0) 



2uj' 



v 2g2 



(17) 



in turns places greatest weight on back-scattering pro- 
cesses with cos(0) = —vq/2oj = —1. (For each q there 
is one relevant fc such that k = — (k — q): in the second 
line we placed the x-axis on the k direction and replaced 
\k\ by u to account for the 5-function.) However the im- 
portance of these processes also depends on the matrix 
element of ero/3 between the initial and final states: 



A 0/3 (|q| = 2lu) oc (* fc |a 0/3 |* fe _ q ) 

= l± e i (<l>k-<t>-h) = \ -J- e i7r 



(18) 



This factor vanishes when the positive sign is taken 
for the Ao case, because pscudospinors with opposite 
momentum in the same valley are orthogonal. The 
same effect is responsible for the Klein-paradox effects 
in graphene transport properties^ and led to a faster 
than usual decay of the Friedel oscillations in real space, 
as calculated by Cheianov and Fal'koi^ and by Bena^. 

As we have just shown, due to pseudospinor-related 
matrix-element effects thc-Ao scattering amplitude across 
a Dirac-cone within the same valley vanishes. In the case 
of A3-scattcring quite the opposite happens and the mo- 
mentum independent terms in the numerator of Eq. (|13|) 
lead to a divergence along \q\ = 2uj. We may therefore 
neglect the non-singular contribution and write: 



As(g) 



dx 
o 

1 

4tt 



dx 



d 2 k -lu 2 + x(l - x)q 2 
(2tt) 2 (fc 2 + A) 2 
-•J 2 + x(l — x)q 2 
ijj 2 + x{l — x)q 2 



, 2uj 



V-z 2 - 1 
and the physical quantity is: 



: arctan ■ 



6Nf(q,w) 



sgn(tj) 
2tt 2 



Im 



V-z 2 - 1 



2iw 



(19) 



(20) 



The above results lead to two very different FT LDOS 
patterns around zero momentum transfer, as shown in 
Fig-©- This may serve as a way to use experimental 
LDOS maps to distinguish between A and A 3 type po- 
tentials in a sample. 



( ) 




FIG. 3: The real and imaginary parts of the FT LDOS pat- 
terns of Ao and A3 type scattering, with momentum transfer 
around K' — K = (— 4n, 0). Features around other inter-valley 
scattering vectors are obtained through 2n/3 rotations. 



B. Inter- valley scattering 

In this subsection we consider the case of momentum 
transfer K 1 — K + q where q is small and therefore the 
initial and final states are each in the vicinity of a valley, 
however not the same valley. This leads to the following 
expression: 



A / 3 (<?) 



d 2 k 



± (k x + ik y )(-k x +q x - ik y + kyL .s 
(w 2 + k 2 ){uj 2 + (q - fc) 2 ) ^ ' 



where the sign in front of k x — q x has been reversed from 
Eq. (ri2|) since the chirality of the final state at fc — q state 
is different from the chirality of the initial state at fc. 
This sign change leads to an angular dependent inten- 
sity peaked along the contour at |g| = 2u. Using the 
Schwinger/Feynman trick with momentum shift as be- 
fore we arrive at: 



A / 3 (<?) 
Ao/3(<?) 



dx 




1 

4tt 



d 2 k — lo 2 ± x(l — x){q x + iq y y 



(2nf 



(k 2 + Af 



2u 



[l±e 2 ^o)G ( g ) ±e 2 ^ 



(22) 



where the function Q(z) is defined as before and <\> q = 
a,rcta,n(q y / q x ) . Note that the physical quantity, 6N(q, ui), 
is obtained by analytic continuation of the frequency and 
has both real and imaginary parts. On the contour de- 
fined by |q| = 2u>, the real part of SN(q) has a sin 2 (0 q ) or 
cos 2 (0q) angular dependence for ao and 03 type scatter- 
ing, respectively. The imaginary part varies as ± sin(20 q ) 
along the contour. These expressions for the two types of 
scattering differ in their orientation, explaining one of the 
principle features of the tight-binding model momentum 
space maps. 

These results are plotted in Fig.©. Note that the so- 
lutions match the ones obtained numerically in the tight 
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FIG. 4: FT LDOS - A zoom in around momentum transfer 
K'-K taken from Fig.{TJ, to be compared with Fig.(|3|- 



binding model (see enlarged features in Fig.Q). In or- 
der to obtain the full momentum space unit cell, however, 
one should apply the phase cxp(ir • Q) to parts of the 
unit cell which are outside of the Brillouin zone. This 
leads to a 2tt/3 rotation about q = 0. 

The LDOS modulations at \q\ = 2lo may be under- 
stood through the dominant scattering contribution as 
before. We may use the same arguments as in the case 
of intra- valley scattering to show that the dominant pro- 



cess in A / 3 (A"' — K + q) has k = — (fe 

id k — q are opposite ex 
valley. The eigenvectors arc now = (l,e~ l ^ h ) 



K), 

so that k and k — q are opposite except for the shift in 

and 

(1, e~ l ^ h -i) where </>' is the angle in the K' val- 
ley. In the K' valley the sign of k x is reversed (and the 
chirality of the bands is reversed) so that the angle be- 
comes <bf. — it — (j>k ■ This leads to 



A 0/3 (q\ = 2w) oc (* fc |a 0/3 |*' fc _ q ) 
= 1 ± e^^^'k-q) = \ ± e 2i ^ q . 



(23) 



The real and imaginary parts of this expression reproduce 
the angular dependence of Fig.©. 



C. Off-Diagonal disorder 

In the previous sections we have modeled an impurity 
as an on-site potential. It is reasonable to assume that 
impurities may induce a change in the hopping ampli- 
tude. Such an effect can be realized, for example, by 
stretching of bonds out of the graphenc plane. In this 
paper we consider a hopping amplitude suppression on 
three bonds around an impurity on the A site. This per- 
turbation is non-local, breaks the sublattice symmetry 
and is described by an off-diagonal potential matrix. 

Let us define the following perturbation Hamiltonian: 



8H 



5>(r)£[ctd r+a +4 +5 c T 



(24) 



where 6t(r) is the change in hopping amplitude around 
the A atom in unit cell r. The operators c and d annihi- 
late an electron on the A and B sublattice respectively. 
Writing the perturbation in momentum space and with 
a pseudospinor vector ^ = {c\ eft) we find: 



T(k,k-q) 



7fc-q 



(25) 



In the Born approximation this leads to G(q,oj) 
Y /k G (k : uj)T(k,k- q)G°(k- q,w) such that 



G AA {q, iuj) = iuj > r 

^ +7fe)(w 2 +lk- q ) 

G B b (q, iui) = 2iuj ^ 



7fc7fc-q 

+ ll){u 2 + lk- q ) 



(26) 



The numerical and analytic evaluation of the above ex- 
pressions is similar to that presented before. We mention, 
however, that this perturbation is a symmetric function 
of the bias voltage to and therefore may be separated 
from the on site potential which produces antisymmet- 
ric functions. Since this perturbation is centered around 
an A atom the LDOS measured on the A sublattices is 
symmetric in real space and does not have an imaginary 
part when Fourier transformed. On the B sublattice, 
the effects of inter-sublattice scattering is seen and one 
is able to see both real and imaginary parts with an an- 
gular dependent amplitude around inter- valley scattering 
vectors. When combining both A and B sublattices the 
result is very similar to Ao(q,(jj) which was presented in 
the previous subsections. In the linear approximation the 
off diagonal contribution is SN ff(q 7 u;) = 2u;5N^-(q,uj). 



V. SUMMARY AND DISCUSSION 

We have obtained analytic expressions for STM 
momentum-space LDOS maps in graphene using the 
massless Dirac equation model for this material. We find 
that smooth disorder produces features near q = and 
reciprocal lattice vectors. The most interesting and sur- 
prising feature is the absence of the back-scattering peak 
which would be expected on the basis of scattering phase 
space considerations. The feature is absent because of 
the sublattice pscudospin chirality of Dirac band states 
which causes disorder induced back-scattering matrix el- 
ements to vanish, the same feature of graphenc which 
helps to enhance its mobility. Atomic length scale dis- 
order leads to both q = features and to features near 
the inter-valley scattering wavevector. For these features, 
pscudospin chirality does not cause backscattering ma- 
trix elements to vanish, and instead leads to an interest- 
ing angular patterning of the on-shcll peak in the LDOS 
map. Our Dirac model analytic results are in agreement 
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with tight-binding model numerical results. These results 
demonstrate the potential of STM experiments to shed 
light on the character of disorder in a graphenc sample. 

In this paper we have not accounted for the influence 
of electron-electron interactions on the tunneling DOS 
of a graphene system. Indeed the potential of LDOS- 
map experiments to provide a high resolution probe of 
interaction effects in the one-particle Greens' function of 
graphenc sheets is a major motivation for undertaking 
these experiments. On the basis of existing theor y 38 i 39 
it appears that electron-electron interactions alter the 
quasiparticle velocity, introduce lifetime broadening, and 
also under some circumstances, introduce sidebands as- 
sociated with plasmon emission. In the present theory 
the LDOS-map depends only on energy relative to the 
Dirac point, whereas in an interacting system the map 
will also depend on the ratio of the energy to the Fermi 
energy. Although the role of interactions must therefore 
be considered carefully in analyzing LDOS-map exper- 
iments, the disorder and matrix element considerations 
explained here will still play a primary role. One role 
played by electron-electron interactions will be that of 



screenin g 27 i the external disorder potential, weakening 
its strength. The weak-disorder theory explained here is 
more likely to be qualitatively correct when the Fermi 
level lies away from the Dirac point, so that the Fermi- 
level density-of-states is larger, the screened potentials 
weake r 27 i 41 , and the Born approximation analysis more 
accurate. The Born approximation analysis is proba- 
bly invalid at small average carrier densities where the 
carrier spatial distribution appear a 42 i 43 highly inhomo- 
geneous. LDOS-map experiments should be interesting 
in both regimes. 



Acknowledgments 

The authors thank C. Bcna for useful comments. This 
work was supported in part by the Welch Foundation, by 
the Army Research Office, by the Natural Sciences and 
Engineering and Research Council of Canada and by the 
NRI-SWAN program. 



1 K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, Y. 
Zhang, S.V. Dubonos, I.V. Grigorieva, A. A. Firsov, Sci- 
ence 306, 666 (2004). 

2 A.K. Geim and K. S. Novoselov, Nat. Mater. 6, 183 (2007). 

3 A.K. Geim and A.H. MacDonald, Physics Today 60 , 35 
(2007). 

4 K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. 
Katsnelson, I.V. Grigorieva, S.V. Dubonos, A. A. Firsov, 
Nature 438, 197 (2005); 

5 Yuanbo Zhang, Yan-Wen Tan, Horst L. Stormer, Philip 
Kim, Nature 438, 201 (2005); Y. Zhang, Z. Jiang, J. P. 
Small, M. S. Purewal, Y.-W. Tan, M. Fazlollahi, J. D. 
Chudow, J. A. Jaszczak, H. L. Stormer, P. Kim Phys. Rev. 
Lett. 96, 136806 (2006). 

6 C.L. Kane and E.J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

7 NA. Sinitsyn, A. H. MacDonald, T. Jungwirth, V.K. 
Dugaev, J. Sinova, Phys. Rev. B 75, 045315 (2007). 

8 P. R. Wallace, Phys. Rev. 71, 622 (1947). 

9 S. H. Pan, J. P. ONeal, R. L. Badzey, C. Chamon, H. 
Ding, J. R. Engelbrecht, Z. Wang, H. Eisaki, S. Uchida, 
A. K. Gupta, K. W. Ng, E. W. Hudson, K. M. Lang, J. C. 
Davis, Nature 413, 282 (2001); C. Howald, P. Fournier and 
A. Kapitulnik, Phys. Rev. B 64, 100504(R) (2001); J. E. 
Hoffman, K. McElroy, D.-H. Lee, K. M Lang, H. Eisaki, S. 
Uchida, J.C. Davis, Science 297 1149 (2002); K. McElroy 
Nature 422, 592 (2003). 

P. Mallet, F. Varchon, C. Naud, L. Magaud, C. Berger, 
and J.-Y. Veuillen Phys. Rev. B 76, 041403(R) (2007); G. 
M. Rutter, J. N. Crain, N. P. Guisinger, T. Li, P. N. First, 
J. A. Stroscio, Science 317 217 (2007); 

1 Cristina Bena and Steven A. Kivelson Phys. Rev. B 72, 
125432 (2005) 

2 Christina Bena, Phys. Rev. Lett. 100, 076601 (2008). 

3 Z. F. Wang, Ruoxi Xiang, Q. W. Shi, Jinlong Yang, Xi- 



aoping Wang, J. G. Hou, and Jie Chen, Phys. Rev. B 74, 
125417 (2006). 

Vadim V. Cheianov and Vladimir I. Falko Phys. Rev. Lett. 

97, 226801 (2006) 

Vadim V. Cheianov, Vladimir Falko and B. L. Altshuler 
Science 315, 1252 (2007). 

N. M. R. Peres, F. D. Klironomos, S.-W. Tsai, J. R. San- 
tos, J. M. B. Lopes dos Santos and A. H. Castro Neto, 
Europhysics Letters 80, 67007 (2007). 
Q.-H. Wang and D.-H. Lee, Phys. Rev. B 67, 020511(R) 
(2003); L. Capriotti, D. J. Scalapino, and R. D. Sedgewick, 
Phys. Rev. B 68, 014508 (2003); L. Zhu, W. A. Atkinson, 
and P. J. Hirschfeld, Phys. Rev. B 69, 060503(R) (2004). 
T. Pereg-Barnea and M. Franz Phys. Rev. B 68, 180506(R) 
(2003). 

T. Pereg-Barnea and M.Franz, International Journal of 
Modern Physics B 19, 731(2005). 

See for example E.N. Economou, Greens Functions in 

Quantum Physics (Springer- Verlag, Berlin, 1979). 

J. Friedel, Nuovo Cimento Suppl. 7, 287 (1958). 

R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical 

Properties of Carbon Nanotubes (Imperial, London, 1998) . 

G.W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 

M. I. Katsnelson, A. K. Geim, Phil. Trans. R. Soc. A 366 

195 (2008). 

Fernando de Juan, Alberto Cortijo and Mara A. H. Voz- 

mediano, Phys. Rev. B 76, 165409 (2007). 

A. Fasolino, J. H. Los, M. I. Katsnels on, larXiv:0704.1793l 
Eun-Ah Kim and A. H. Castro Neto, |cond-mat /0702562] 
Kentaro Nomura, A.H. MacDonald, Phys. Rev. Lett. 96, 
256602 (2006). 

E.H. Hwang, S. Adam and S. Das Sarma, Phys. Rev. Lett. 

98, 186806 (2007). 

Vitor M. Pereira, Johan Nilsson, A. H. Castro Neto, Phys. 
Rev. Lett. 99, 166802 (2007). 



9 



30 Vitor M. Pereira, J. M. B. Lopes dos Santos, A. H. Castro 
Neto, larXiv:0712.0806l 

31 T. O. Wehling, A. V. Balatsky, M. I. Katsnelson, A. I. 
Lichtenstein, K. Scharnberg, and R. Wiesendanger Phys. 
Rev. B 75, 125425 (2007). 

32 Cristina Bena and Gilles Montambaux, arXiv:0712.0765 

33 Note that although the the off-diagonal Greens' functions 
are gauge dependent, G\]\r, r,ui) is not. It is important in 
evaluating (r , r , w) to be careful that the gauge choices 
made in evaluating G' ' are applied consistently, particu- 
lar the choice made for the relative phase of Bloch state 
wavefunction amplitudes on different sublattices. 

34 Yafis Barlas, T. Pereg-Barnea, Marco Polini, Reza Asgari 
and A.H. MacDonald Phys. Rev. Lett. 98 236601 (2007). 

35 Marco Polini, Reza Asgari, Yafis Barlas, T. Pereg- 
Barnea and A.H. MacDonald Solid State Communications 
2007.04.035. 

36 See e.g. M. E. Peskin and D. V. Schroeder, An Introduction 



to Quantum Field Theory, Addison- Wesley 1995. 

M. I. Katsnelson, K. S. Novoselov and A. K. Geim, Nature 

Phys. 2, 620-625 (2006). 

Marco Polini, Reza Asgari, Giovanni Borghi, Yafis Barlas, 
T. Pereg-Barnea, and A. H. MacDonald Phys. Rev. B 77, 
08141 1(R) (2008). 

E.H. Hwang and S. DasSarma, Phys. Rev. B 77, 081412(R) 
(2008). 

M. M. Fogler, D. S. Novikov, and B. I. Shklovskii, Phys. 
Rev. B 76, 233402 (2007). 

Y.-W. Tan, Y. Zhang, K. Bolotin, Y. Zhao, S. Adam, E. H. 
Hwang, S. Das Sarma, H. L. Stormer, and P. Kim, Phys. 
Rev. Lett. 99, 246803 (2007);J. H. Chen, C. Jang, M . S. 
Fuhrer, E. D. Williams, M. Ishigami, larXiv:0708.2408l 
J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, J. H. 
Smet, K. von Klitzing and A. Yacoby arXiv:0705.2180 
Victor M. Galitski, Shaffique Adam and S. Das Sarma, 
Phys. Rev. B 76, 245405 (2007). 



